Community composition, bacterial symbionts, antibacterial and antioxidant activities of honeybee-associated fungi

Background Fungi associated with insects represent one potentially rich source for the discovery of novel metabolites. However, a comprehensive understanding of the fungal communities of Apis mellifera ligustica remains elusive. Results Here, we investigated the phylogenetic diversity and community composition of honeybee-associated fungi using combination of culture-dependent and culture-independent approaches. A total of forty-five fungi were isolated and purified from the Apis mellifera ligustica, royal jelly, and honeycomb, which belonged to four classes and eleven different genera. Furthermore, 28 bacterial 16S rRNA gene sequences were obtained by PCR from the fungal metagenome. High-throughput sequencing analyses revealed that the fungal communities were more diverse, a total of 62 fungal genera were detected in the honeybee gut by culture-independent method, whereas only 4 genera were isolated by culture-dependent method. Similarly, 247 fungal genera were detected in the honeycomb, whereas only 4 genera were isolated. In addition, we assessed the antibacterial and antioxidant activities of fungal isolates. Most fungal crude extracts obtained from the cultivation supernatant exhibited antioxidant activities. Only two fungal crude extracts displayed moderate activity against Escherichia coli and Staphylococcus aureus. Chemical analysis of Chaetomium subaffine MFFC22 led to the discovery of three known compounds, including cochliodinol (1), emodin (2), chrysophanol (3). Among them, cochliodinol (1) showed intense DPPH radical scavenging activity with the 50% inhibitory concentration (IC50) of 3.06 μg/mL, which was comparable to that of the positive ascorbic acid (IC50 = 2.25 μg/mL). Compound 2 displayed weak inhibitory activities against Micrococcus tetragenus and S. aureus. Conclusions This research provided a fundamental clue for the complex interactions among honeybees, fungi, bacterial symbionts, and the effects on the honeybee. Furthermore, the diversity of honeybee-associated fungi had great potential in finding the resource of new species and antioxidants. Supplementary Information The online version contains supplementary material available at 10.1186/s12866-022-02580-4.

and bees [3][4][5]. It is worth mentioning that the relationship between insects and their associated fungi capture researchers' attention, for instance, beetle-fungus farming symbiosis, fungus-cultivating termites, and fungusfarming ants [6][7][8]. In addition, reports indicate that some bacteria were present in fungal hyphae [9,10]. It is worth considering to evaluate in many further aspects of the associations between fungi and insects as well as fungi and bacteria.
Insect-associated fungi are vital microbial sources of natural bioactive products. The insect-associated fungi isolated from arthropod cadavers, leaf-cutting ants, and stoneflies were contained new microbial species [11][12][13]. Insect-associated fungi could enhance insects' fitness by producing bioactive compounds [14]. Furthermore, many novel compounds have been discovered from insect-associated fungi, and these compounds had the potential as immune inhibitors, antibacterial agents, and biofungicides [15][16][17]. However, the current research in metabolites produced by insect-associated fungi is still not comprehensive.
As a social insect, honeybees belong to the class Insecta and the order Hymenoptera in the phylum Arthropoda. It is a critical species for agricultural production as pollinators [18]. Various parts of honeybee larvae and adults, their food, and honeycomb harbor numerous microorganisms, which play a significant role in food digestion, pollination, and antagonistic effect against different pathogens [19]. Especially, fungi associated with honeybee may provide material for pollen degradation or assist in royal jelly maturation, and also can be a food source [20,21]. For example, the Brazilian stingless bee larva grows by eating fungal mycelia of Monascus inside the brood cell [22]. Therefore, honeybees are a potential model of fungus-host-symbiont interactions, which is worth exploring. As a high royal jelly-producing honeybee, Apis mellifera ligustica acquired fungi from their diet, the surrounding environment, or mates [23]. However, only a few reports involve fungi isolated from specimens associated with A. mellifera ligustica [24]. Our knowledge about the species, biological activity, and secondary metabolites of these associated fungi are limited until now. Here, we investigated the diversity of fungi from A. mellifera ligustica using culture-dependent and cultureindependent methods and explored the bacterial symbionts, biological activity, and secondary metabolites of the fungi isolated from A. mellifera ligustica (larvae, adults), honeycomb, and royal jelly.

Identification of cultivable fungi associated with honeybee
In this study, a total of forty-five fungi were isolated from the honeycomb, royal jelly, larvae, and different parts of honeybee on nine different media (Table 1, Fig. 1). The morphological and microscopic characteristics of the ten representative strains are shown in Fig. 2. Among them, seven strains were isolated from honeycomb, thirteen from royal jelly, ten from larvae, two from honeybee cuticle, seven from gut, three from head, and three from hypopharyngeal gland. ITS sequence analysis of the forty-five isolates revealed that the fungal isolates belonged to two different phyla (Basidiomycota and Ascomycota) and eleven different genera. Twenty-seven isolates (60.0%) were distributed in the Agaricomycetes within the phylum Basidiomycota. The other eighteen isolates were grouped into three classes [Sordariomycetes (26.7%), Dothideomycetes (11.1%), and Eurotiomycetes (2.2%)] within Ascomycota.
Fungal isolates belonging to class Agaricomycetes were assigned to two orders, including the Agaricales (3 isolates) and Polyporales (24 isolates). The three strains in Agaricales were isolated from the honeycomb, royal jelly, and hypopharyngeal gland, respectively, which showed highly similar to Schizophyllum commune with more than 99% identity. The three strains in Agaricales were isolated from the honeycomb, royal jelly, and hypopharyngeal gland, respectively, which showed highly similar to Schizophyllum commune with an identity of more than 99%. The strains of Polyporales were isolated from most of the samples and represented by the family Polyporaceae, Irpicaceae and Phanerochaetaceae. Among them, most strains (17 isolates) belonging to genus Irpex were identified as Irpex lacteus. Notably, the strain FWJ13 showed only 98% similarity to I. lacteus. In addition, two strains isolated from honeybee larvae were identified as Bjerkandera adusta and Leiotrametes lactinea, respectively. Moreover, B. adusta was also found in royal jelly. The last four strains belonging to Trametes genus showed similar to Trametes hirsuta with an identity of 99%.
The other representative class was Sordariomycetes, including three orders (Hypocreales, Sordariales, and Xylariales). Eight isolates belonging to the family Hyporcreaceae exhibited a sequence match of more than 99% to Fusarium solani. Fusarium spp. were found in association with honeybee gut, larvae, and royal jelly. Both strains belonging to Sordariales were isolated from the honeycomb, identified as Chaetomium subaffine with an identity of 100%. The other two strains belonging to the Xylariales were identified as Arthrinium arundinis, which were found to be most related to honeybee intestines and larvae.
The five fungal sequences of Dothideomycetes were grouped into Pleosporales, including Paraconiothyrium brasiliense (3 isolates), and Alternaria alternate (2 isolates). The P. brasiliense with an identity of 99% were isolated from honeybee cuticle and larvae, respectively. The genus Alternaria was only isolated from honeybee larvae. Finally, only one isolate belonging to class Eurotiomycetes was grouped into the order Eurotiales, which belonged to the genus Aspergillus with a high sequence match to Aspergillus flavus (> 99%). The strain was isolated from honeybee gut.

Culture-independent community
The ITS1 region was performed to analyze the fungal community within the honeybee gut and honeycomb by using Illumina Miseq sequencing. A total of 408,112 high-quality fungal clean reads were generated from two populations with 3 replicates for each population (honeybee gut: MFF1_1, MFF1_2, MFF1_3, honeycomb: FC_01, FC_02, FC_03) (Supplementary Table S1 and S2). A total of 3 fungal phyla were identified in honeybee gut samples, including Ascomycota, Basidiomycota, and unclassified fungi, the average abundance was 93.29, 6.67, 0.04%, respectively (Fig. 3A). Ascomycota (88.81% of the average abundance) and Basidiomycota (7.05%) were identified in honeycomb samples (Fig. 3C). Regarding the composition of fungal community, all samples were dominated at the phylum level by Ascomycota. The proportions of Ascomycota in the honeybee gut were higher than in the honeycomb. The community composition was further analyzed at the genus level. 62, 247 genera were identified across investigated samples of honeybee gut and honeycomb, respectively. Among them, Kodamaea (80.13%), Zygosaccharomyces (8.19%), Wallemia (6.11%), and Wickerhamomyces (3.62%) had higher abundance in the honeybee gut (Fig. 3B). However, Bipolaris (47.07%), Metschnikowia (18.17%), Starmerella (3.85%), Trichoderma (16.77%), Kodamaea (3.35%) had higher abundance in the honeycomb (Fig. 3D). The relative abundance of yeast such as Kodamaea, Zygosaccharomyces, and Wickerhamomyces was higher in the honeybee gut than in the honeycomb. There are differences in fungal communities of honeybee gut and honeycomb samples between the two approaches. Culture-independent community analysis showed more diverse fungal communities than did the culture-dependent method.

Identification of bacterial symbionts
In this study, a total of twenty-eight bacterial symbionts 16S rRNA were observed from genomic DNA extractions of the honeybee-associated fungi (Table2). Among them, nineteen, six, two, and one strain were isolated from the host of the Agaricomycetes, Sordariomycetes, Dothideomycetes, and Eurotiomycetes, respectively. All bacterial symbionts were divided into four genera by initial BLAST comparisons in GenBank and were putatively identified as Bacillus velezensis, Bacillus siamensis, Microbacterium tumbae, Pandoraea sputorum, and Achromobacter xylosoxidans. These bacterial symbionts were mainly identified as B. velezensis and B. siamensis in Bacillaceae. Moreover, B. velezensis was obtained from fungi of different honeybee samples or parts, such as honeybee gut, hypopharyngeal gland, head, cuticle, larvae, honeycomb, and royal jelly.

Antimicrobial activities
The filter paper disc method was used to evaluate the antibacterial activities of 36 fungal crude extracts, which were obtained from the cultivation supernatant (Table3). Most of the fungal crude extracts showed no antibacterial activities. However, I. lacteus FWJ16 displayed moderate activity against E. coli with a disc diameter of inhibition zone diameter (IZD) of 10.13 mm, which was weaker than the positive gentamicin sulfate with the IZD of 26.33 mm. Additionally, P. brasiliense MFYC10 also displayed moderate activity against S. aureus with the IZD of 9.67 mm, which was weaker than the positive gentamicin sulfate with the IZD of 23.92 mm.

Antioxidant activities
The antioxidant activities of 36 fungal crude extracts are shown in Table3. The results revealed that 32 extracts (88.9%) exhibited antioxidant activities under the concentration of 166.67 µg/mL. Among them, scavenging rates between 10-40%, 40-70%, and 70-100% were found in 19 (52.8%), 3 (8.3%), and 5 (13.9%) strains, respectively. Especially, crude extract of MFFC22 showed the strongest antioxidant activity with the DPPH scavenging activity of 91.6%. Thus, the MFFC22 was further selected to  Fig.4). The result showed that the MFFC22 had a higher activity with the increasing concentration. Its IC 50 on DPPH scavenging activity was 22.11 µg/mL, which was relatively weaker than a commercial natural antioxidant ascorbic acid (Vc, 2.25 µg/mL).

Antibacterial and antioxidant activities of compounds
The antibacterial activities of the compounds isolated from MFFC22 were shown in Table 4. Specifically, compound 2 displayed weak inhibitory activities against M. tetragenus and S. aureus with the IZD of 7.00 and 7.33 mm, respectively. Compound 1 had no inhibition effect on three bacterial strains.
The antioxidant activities of compound 1 isolated from MFFC22 are shown in Fig. 4. Notably, the 50% inhibitory concentration (IC 50 ) of cochliodinol (1) on DPPH scavenging activity was 3.06 μg/mL, which was comparable to that of positive ascorbic acid (2.25 μg/mL).

Discussion
Microbial community is extensive among insects. As a special kind of microbial community, insect-associated fungi have been a source of new microbial resources [29,30] and biologically active natural products [16,31,32]. While fungi were present in A. mellifera ligustica, the leading research mainly focused on the honeycomb, honeybee gut, pollen, and bread [23,[33][34][35][36]. To the best of our knowledge, this is the first report on the diversity, antibacterial and antioxidant activities of culturable fungi from honeycomb, honeybee product (royal jelly), larvae, and different parts of A. mellifera ligustica. Forty-five fungi were isolated and characterized by the dilution-plate method and molecular biological identification. Furthermore, three metabolites were purified and characterized from C. subaffine MFFC22. Consequently, honeybee-associated fungi can provide a resource for microbial diversity and natural products.
The same fungus of I. lacteus was isolated and identified from honeybee gut, head, hypopharyngeal gland, larvae, and royal jelly in our study. It suggests that the fungus is widely distributed in A. mellifera ligustica. Some fungi isolated from honeybee larvae, hypopharyngeal gland, honeycomb, and royal jelly were wood-rotting fungi, such as T. hirsuta, S. commune, B. adusta, and L. lactinea [37][38][39]. Besides, the genus Arthrinium isolated from honeybee gut and larvae has been reported as a plant pathogen [40]. Compared with the fungal diversity of other insects, some fungi isolated from honeybee were the same as those of other different insects. For example, I. lacteus and P. brasiliense were also isolated from termites, and F. solani was isolated from beetles [41,42]. Aspergillus, Alternaria, Chaetomium, Fusarium, and Arthrinium genera have already been found in the honeybee. For example, Aspergillus, Chaetomium, Fusarium, and Alternaria sp. were isolated from honeybee gut, and Arthrinium sp. was isolated from pollen [23,[33][34][35]. However, most of these fungi reported in this study are new for this honeybee. For instance, this is the first report of the isolation of T. hirsuta and S. commune from both honeycomb and royal jelly, I. lacteus from all parts of honeybee samples. Additionally, one rare fungus P. brasiliense was found in honeybee cuticle, which was reported previously in the gut of Acrida cinerea [43].
A total of 62 fungal genera were discovered in honeybee gut by culture-independent method. Among them, Kodamaea, Zygosaccharomyces, Wallemia, and Wickerhamomyces were predominant genera. Previous study  [20]. However, Starmerella and Hanseniaspora had high relative abundance in honeybee gut samples collected from Italy and Saudi Arabia [16]. Previous studies on predominant genera were inconsistent with each other, and also differed from our results. Presumably, it was due to the different surrounding environments of the sampling sites. Some genera analyzed by culture-independent approaches were also isolated by culture-dependent method, such as Aspergillus in honeybee gut, Trametes, and Chaetomium in honeycomb. However, cultureindependent approaches in general revealed higher microbial composition and diversity compared to culture-dependent method. For example, 247 fungal genera were detected in the honeycomb by culture-independent method, whereas only 4 genera were isolated. Interestingly, as the causative agent of chalkbrood disease in honeybee larvae [44], Ascosphaera spp. was detected in honeybee gut and honeycomb by culture-independent approach, but it has not been isolated yet. The DPPH radical scavenging assay indicated that most crude extracts from honeybee-associated fungi showed antioxidant activities. Notably, the crude extract of C. subaffine MFFC22 isolated from honeycomb has shown notable antioxidant activity. Chaetomium genus is already known to have strong DPPH scavenging activity, such as C. nigricolor, C. globosum, C. cruentum [26,44,45]. Nevertheless, C. subaffine and its metabolites have not reported DPPH scavenging activity. The compound 1 from MFFC22 showed strong DPPH scavenging activity with the IC 50 values of 3.06 μg/mL, comparable to the positive ascorbic acid (Vc). Comparing a commercial antioxidant BHT with IC 50 of 95.7 μM, Dehghan et al. found that compounds 2 and 3 also exerted moderate antioxidant activities against DPPH with IC 50 of 271.2 and 297.0 μM, respectively [46]. The phenolic compounds with the free hydroxyl groups were well known for robust antioxidant activity [47]. Therefore, the antioxidant activity of compounds 1-3 may be attributed to their phenolic hydroxyl groups. Additionally, it was reported that antioxidant compounds produced by fungi could protect their hosts by enhancing tolerance to abiotic stresses [48]. Whether compounds 1-3 have similar effects deserves further exploration. In addition, these isolated compounds have been shown to have other bioactivities. For example, compound 1 was a common fungal metabolite with cytotoxic activities against the KB, MDA-MB-435 and MRC5 cell lines [49]. Compounds 2-3 were anthraquinone derivative and acted as an antimalarial and antiallergic agents [50,51].
Symbiotic bacteria have occurred in some fungi, which affected fungi function and subsequent host-fungus interactions. However, only a few studies have explored their effects in insect-associated fungi [52,53]. Here, we found the first insight into the bacterial symbionts in fungi associated with honeybee. The bacterial symbiont  Gentamicin sulfate a 26.00 ± 0.00 20.00 ± 0.00 19.00 ± 0.00 from honeybee has already been shown to protect the host by producing inhibitory compounds or providing nutrition [54,55]. However, whether fungal-bacterial symbionts influence the growth of honeybee-associated fungi, and how they affect tripartite honeybee-fungusbacterium mutualisms, remains to be explored further.

Conclusions
Here, the diversity of the honeybee-associated fungi was investigated using both culture-dependent and cultureindependent analysis methods. This study expands our knowledge of honeybee-associated fungi and further raises the pool of fungal species from A. mellifera ligustica. The results show that several of these fungi have antibacterial and antioxidant activities, among which the fungus C. subaffine MFFC22 was the most prominent. Furthermore, the antioxidant activities of C. subaffine MFFC22 can be attributed to the identified phenolic compounds. Collectively, the culturable honeybee-associated fungi provide insight into the widespread insect symbionts, which have great potential in finding the resource of novel species and antioxidants.

Sample collection and microbial isolations
A. mellifera ligustica (including larvae, adults), honeycomb, and royal jelly were collected from the Institute of Apicultural Research, Anhui Agricultural University, Hefei, China (GPS: 31 • 53ʹN, 117 • 20ʹE) between April and July 2021. The honeybee larvae and adults starved for 24 h, and all samples were stored at 4 ℃. Initially, seven larvae, seven adults, and one gram honeycomb were placed separately into 10 mL sterile phosphate-buffered saline (PBS) buffer solution in an autoclaved 50 mL tube to obtain fungi from external isolation. Then, the same tube from the external was filled with 10 mL 75% ethanol for 2 min, followed by rinsing in sterilized 10 mL 1% bleach with 0.1% tween 20 for three times (30 s each). The supernatant was removed and replaced with 10 mL sterile PBS solution. Subsequently, sterile forceps were used to dissect samples of adult honeybee to get head, gut, and hypopharyngeal gland. Each body part was grounded in 10 mL sterile PBS. According to the earlier report, the honeybee larvae, one gram honeycomb, and one gram royal jelly were fully homogenized separately in 10 mL sterile PBS [56]. Then, the homogenates were diluted in a tenfold series (i.e., 10 −1 , 10 −2 , 10 −3 ), and aliquots of 100 µL from each dilution were spread onto nine isolation media (Table 5). Pure colonies of fungi from the appropriate dilution were transferred into a new PDA medium and incubated aerobically at 28 ℃. All these isolated fungal strains were preserved on PDA slants at 4 ℃ until use. The fungi were used for freeze-drying preservation by a freeze-dryer (BTP-3ES; SP Scientific, USA) [57], and were stored at our institute.

DNA sequencing
All honeybee-associated fungi were identified by molecular techniques and morphological characteristics [10,58,59]. Each fungus was cultured in malt extract medium at 28 ± 0.5 ℃ for 7 days. Then the fungal genomic DNA was extracted using the Fast DNA Extraction Kit (Aidlab Biotechnologies Co., Ltd., Beijing, China) as claimed by the manufacturer's specification. The primers ITS1/ITS4 were used to amplify ITS based on the fungal genomic DNA. Additionally, 16S rRNA was amplified using the primers 27F/1492R by PCR from the fungal metagenome [10,60]. The quality of PCR products was visualized on 1% agarose gel by electrophoresis. Each product Table 5 Culture media for the isolation of honeybee-associated fungi was successfully amplified from a PCR for sequencing (Tsingke Biotechnology Co., Ltd., Beijing, China).

Identification of fungi and bacterial symbionts
As mentioned before [42], all achieved sequences' affiliation returned from Tsingke Biotechnology Company recognized by the available data in BLAST from the National Center for Biotechnology Information (NCBI) database. Sequence alignment and Neighbor-joining phylogenetic analysis were performed using MEGA software version 5.0. Bootstrap analysis of tree construction on the strength of the sequences was habituated to judge the neighbor-joining information based on 1,000 replicates [61]. The obtained ITS sequences were deposited in the GenBank database under accession numbers OK184563-OK184606 and OK285068. The accession numbers of bacterial 16S rRNA sequences were OK147622-OK147645 and OK169608-OK169611.

Culture-independent community analysis
The sample pretreatment was the same as those mentioned above to obtain the honeybee gut and honeycomb. Then, the total genome DNA of samples was extracted using the Fast DNA Extraction Kit, the concentration and purity of DNA were confirmed on 2% agarose gels. Each sample was performed in triplicate. ITS1 genes of regions were amplified using a specific primer (ITS5-1737F and ITS2-2043R) with the barcode. The reaction conditions were 98℃ for 1 min, followed by 30 cycles of 98℃ for 10 s, 50℃ for 30 s, and 72℃ for 30 s, finally 72℃ for 5 min. Mixture PCR products were purified using Qiagen Gel Extraction Kit (Qiagen, Germany). Sequencing library was generated using TruSeq ® DNA PCR-Free Sample Preparation Kit (Illumina, USA) according to manufacturer's recommendations and index codes were added. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100 system. Finally, the library was sequenced on an Illumina NovaSeq platform using 250 bp pairedend reads. Raw reads were demultiplexed and quality-filtered using QIIME V1.9.1 [62]. Effective Tags were obtained by comparing with the reference Unite database and using UCHIME Algorithm to remove chimera sequences [63,64]. Sequences with ≥ 97% similarity were assigned to the same operational taxonomic units (OTUs) by using the software Uparse V7.0.1001 [65]. Taxonomy of each representative sequence was assessed by using the MUS-CLE software V3.8.31, and comparison against the Unite database based on blast algorithm [66]. Raw data is available from the NCBI Short Read Archive under accession number PRJNA817087 and PRJNA817099.

Cultivation of fungi and preparation of culture extract
Each fungus was inoculated in PDA medium and incubated at 28 ± 0.5 ℃for 3-4 days. Then, the fresh mycelia of each fungus (9 plugs of 5 mm) were inoculated in a 250 mL Erlenmeyer flask containing 150 mL of ME liquid medium and incubated at 28 ± 0.5 ℃ for 7 days in a shaker rotating at 180 rpm. The culture was passed through four layers of cheesecloth to remove the fungal thallus, then the supernatant was exhaustively extracted three times with ethyl acetate (EtOAc, 1:1, v/v). The fungal crude extracts were concentrated in a vacuum to yield culture extract for further experimental use.

Antimicrobial activities
The filter paper disc method was used to screen the antibacterial activity of crude extracts [67]. Three bacterial strains (Escherichia coli (ATCC8739), Micrococcus tetragenus (ATCC35098), and Staphylococcus aureus (ATCC6538) were selected for the test and cultivated on tryptic soy blood agar (TSBA) medium in 37 ℃. Next, sterile filter paper discs (5 mm in diameter) were added to 5 µL of the tested crude extracts, which were dissolved separately in acetone to get a concentration of 18 mg/mL. The filter paper disc treated with acetone alone and gentamicin sulfate was served as a negative and positive control, respectively. The plates were prepared in triplicate and were cultivated at 37 ℃ for 24-36 h. Finally, the diameters of inhibition zone diameter (in mm) were measured for evaluating antimicrobial activity.

Antioxidant activities
The 2,2-diphenyl-1 picrylhydrazyl (DPPH) radical scavenging activity was conducted according to the previous method with some modifications [68]. 1 mL of each crude extract (1 mg/mL) mixed with 5 mL of DPPH solutions in methanol (20 μg/mL). The mixture was incubated at room temperature for 30 min in darkness, and the absorbance was measured at 517 nm using a spectrophotometer (UV-1601; Beijing Beifenruili Analytical Instrument Co., Ltd., China). Methanol and ascorbic acid (Vc) were used as negative and positive controls, respectively. Each sample was performed in triplicate. The radical scavenging activity of crude extracts in the DPPH was calculated as follows: Scavenging rate (SR) (%) = (A b -A s ) /A b *100% where SR was the scavenging activity of the tested sample (%), A b was the absorbance without sample, and A s was the absorbance in the presence of the samples or a positive substance.

Structural elucidation of metabolites
The structures of all compounds were primarily analyzed by High Resolution-Mass Spectrometry (HR-ESI-MS) and 1 H/ 13 C-nuclear magnetic resonance (NMR) spectroscopies. HR-ESI-MS spectra were measured on a TripeTOF 4600 mass analyzer (Bruker, USA). 1 H/ 13 C NMR data were acquired using Agilent DD2 600 Hz spectrometer (Agilent, USA) with tetramethylsilane (TMS) as an internal standard, and chemical shifts (δ) were reported as parts per million (ppm) values.

Antibacterial and antioxidant activities of compounds
The raw data supporting the conclusions of this manuscript will be made available to the authors, without undue reservation, to any qualified researcher. Publicly available datasets were analyzed in this study. This data can be found here: The obtained ITS gene sequences were deposited in the GenBank database under accession numbers OK184563-OK184606 and OK285068. The obtained 16S rRNA gene sequences were deposited in the GenBank database under accession numbers OK147622-OK147645 and OK169608-OK169611. Data of high-throughput sequencing by paired-end Illumina technology of ITS 1 gene amplicons can be retrieved from the NCBI Short Read Archive under accession number PRJNA817087 and PRJNA817099.